Curvature induced optical phonon frequency shift in metallic carbon nanotubes 
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The quantum corrections to the frequencies of the F point longitudinal optical (LO) and transverse 
optical (TO) phonon modes in carbon nanotubes are investigated theoretically. The frequency shift 
and broadening of the TO phonon mode strongly depend on the curvature effect due to a special 
electron-phonon coupling in carbon nanotubes, which is shown by the Fermi energy dependence of 
the frequency shift for different nanotube chiralities. It is also shown that the TO mode near the 
r point decouples from electrons due to local gauge symmetry and that a phonon mixing between 
LO and TO modes is absent due to time-reversal symmetry. Some comparison between theory and 
experiment is presented. 



I. INTRODUCTION 

In the Raman spectra of a single wall carbon nanotube 
(SWNT), the two in-plane optical phonon modes, that 
is, the longitudinal optic (LO) and transverse optic (TO) 
phonon modes at the T point in the two-dimensional Bril- 
louin zone (2D BZ), which are degenerate in graphite and 
graphene, split into two peaks, G + and G~ peaks, respec- 
tively^-^ The splitting of the two peaks for SWNTs is 
inversely proportional to the square of the diameter dt 
of SWNTs due to the curvature effect, in which G + does 
not change with changing dt, but the G~ frequency de- 
creases with decreasing dt^^ In particular, for metallic 
SWNTs, the G~ peaks appear at a lower frequency than 
the G~ peaks for semiconducting SWNTs with a similar 
diameter £ The spectra of G~ for metallic SWNTs show 
a much larger spectral width than that for semiconduct- 
ing SWNTs. Further, the spectral G~ feature shows an 
asymmetric lineshape as a function of frequency which is 
known as the Breit- Wiener- Fano (BWF) lineshape,! The 
origin of the BWF lineshape is considered to be due to 
the interaction of discrete phonon states with continuous 
free electron states. 

It has been widely accepted that the frequency shift of 
the G-band is produced by the electron-phonon (el-ph) 
interaction ! 8 ' 9 ' 10 ! 11 An optical phonon changes into an 
electron-hole pair as an intermediate state by the el-ph 
interaction. This process gives the phonon a self-energy. 
The phonon self-energy is sensitive to the Fermi energy, 
Ep. In the case of graphite intercalation compounds in 
which the charge transfer of an electron from a dopant 
to the graphite layer can be controlled by the doping 
atom and its concentration, Eklund et al. observed a 
shift of the G-band frequency with an increase of the 
spectral widths In this case the frequency shifted spec- 
tra show that not only the LO mode but also the TO 
mode are shifted in the same fashion by a dopant. For a 
graphene monolayer, Lazzeri et al. calculated the Ep de- 
pendence of the shift of the G-band frequency? 9 - The LO 
mode softening in metallic SWNTs was shown by Dubay 



et al.) 13 ' 14 on the basis of density functional theory. Re- 
cently Nguyen et al^ and Farhat et al*ZL observed the 
phonon softening effect of SWNTs as a function of Ep by 
field effect doping and electro-chemical doping, respec- 
tively, and their results clearly show that the LO phonon 
modes become soft as a function of Ep. Ando discussed 
the phonon softening for metallic SWNTs as a function 
of the Ep position, in which the phonon softening occurs 
for the LO phonon mode and for a special range of Ep, 
that is, for \E F \ < h^oj^XL 

In this paper, we show that the T point TO phonon 
mode becomes hard when \Ep\ < 7iu;to/2 and has a con- 
siderable broadening for metallic zigzag nanotubes. The 
occurrence of the phonon hardening for the TO mode 
is due to the curvature effect, a special character of 
the electron-phonon coupling, and a basic consequence 
of second-order perturbation theory. We show using a 
gauge symmetry argument that the electrons completely 
decouple from the TO mode near the T point. Besides, 
we show that for a chiral nanotube, both the LO and TO 
modes are softened due to the fact that the direction of 
the TO phonon vibration is not parallel to the nanotube 
circumferential direction^- 8 - Another interest of ours is 
the mixing of LO and TO phonons to form degenerate 
phonon frequencies. When the LO phonon mode be- 
comes soft, a crossing of the LO mode with the TO mode 
occurs at a certain value of Ep. For such a mode crossing, 
we should generally consider the electron-phonon cou- 
pling for degenerate phonon modes to sense the crossing 
or anti-crossing of the two phonon frequencies as a func- 
tion of Ep. We will show by an analytical calculation 
that there is no mixing between the LO and TO phonon 
modes for any case due to time reversal symmetry. 

The organization of the paper is as follows. In Sec. |TT] 
we show our method of calculation and present the re- 
sults for armchair and metallic zigzag SWNTs. In Sec. lIIII 
using effective mass theory, we show how the el-ph inter- 
action depends on the chiral angles of SWNTs, and in 
Sec.|lV]a discussion based on gauge symmetry and time 
reversal symmetry for the el-ph coupling is given. In 
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Sec. El a comparison with the experiments and summary 
are given. 



II. PHONON FREQUENCY SHIFT 

The frequency shift of the L-point LO and TO phonon 
modes for metallic SWNTs is calculated by second-order 
perturbation theory. The phonon energy including the 
cl-ph interaction becomes huj\ = fuv x °^ + hujf 1 (A = 
LO, TO) where u> x is the original phonon frequency 

(2) 

without the el-ph interaction and hto x is given by 

n J2) _ 2 y |(eh(k)|H in tK)| 2 

X V &4 0) - (E e (k) - E h (k) ) + IT A 

x (f(E h (k) - E F ) - f(E e (k) - E F )) . (1) 

The factor 2 in Eq. ([T]) comes from spin degeneracy. hui x ' 
is the quantum correction to the phonon energy due to 
the electron- hole pair creation as shown in Fig. [TJa). In 
Eq. ([1]), (eh(k)|7ii llt |wa) is the matrix element for creat- 
ing an electron-hole pair at momentum k by the el-ph 
interaction, 7ii n t, E e (k) (Eh(k)) is the electron (hole) 
energy and T\ is the decay width. In Fig. [lja), an in- 
termediate electron-hole pair state that has the energy 
of E = E e (k) - E h (k) is shown. We need to sum (X)k) 
over all possible intermediate electron-hole pair states in 
Eq. ([T]), which can have a much larger energy than the 

phonon (E ^> huf^). 

Since (eh(k)|7ii nt |o;A) is a smooth function of E = 
E e (k) — Eh(k) in the denominator of Eq. ([1]), the con- 

(2) 

tribution to hto x from an electron-hole pair depends 
on its energy. In Fig. []Jb), we pl°t the real part and 
imaginary part of the denominator of Eq. Jl}, h(E) = 
l/(hu)(°) — E + iT) as a function of E in the case of 
huj^ = 0.2 eV and r = 5 meV. Here Re{h{E)) is a 
positive (negative) value when E < fwj^ ' (E > Tiuj^) 
and the lower (higher) energy electron-hole pair makes 

(2) 

a positive (negative) contribution to hiv x . Moreover, 
an electron-hole pair satisfying E < 2\E F \ can not con- 
tribute to the energy shift (shaded region in Fig. QJa) 
and (b)) because of the Fermi distribution function f(E) 
in Eq. ([1]). Thus, the quantum correction to the phonon 
energy by an intermediate electron-hole pair can be con- 
trolled by changing the Fermi energy, E F . For example, 
when \E F \ = then huj x 2) takes a minimum value 

at zero temperature since all positive contributions to 

(2) 

Tiuj x are suppressed in Eq. ((T|). Since Re(h(E)) « —l/E 
for E > hJ°\ all high energy intermediate states con- 
tribute to phonon softening if we include all the electronic 
states in the system. Here we introduce a cut-off energy 
at E c = 0.5 eV as Y^ e ^ <Ec in order to avoid such a 
large energy shift in Eq. ([T]). The energy shift due to 
the high-energy intermediate states (Ylu" >E °) can t> e 
neglected by renormalizing hcu^ so as to reproduce the 



experimental results of Raman spectral since the con- 
tribution from £e(k) > E c just gives a constant energy 
shift to hu)( 2 \ We have checked that the present results 
do not depend on the selection of the cut-off energy since 
E c is much larger than hu!^°\ 

Im(h(E)) is nonzero only very close to E — hu>^°\ 
which shows that the phonon can resonantly decay into 
an electron-hole pair with the same energy. It is noted 
that when \E F \ > hu>^°' /2, ~ at zero temperature 
while T\ may take a finite value at a finite temperature. 
In this paper, we calculate T\ self-consistently by calcu- 



lating T A = -lm(TiuJ x 2) ) in Eq. |T|). 

(a) (b) 




2\E F \ E (eV) 

FIG. 1: (color online) (a) An intermediate electron-hole pair 
state that contributes to the energy shift of the optical phonon 
modes is depicted. A phonon mode is denoted by a zigzag 
line and an electron-hole pair is represented by a loop. The 
low energy electron-hole pair satisfying < E < 2\E F \ is 
forbidden at zero temperature by the Pauli principle, (b) The 
energy correction to the phonon energy by an intermediate 
electron-hole pair state, especially the sign of the correction, 
depends on the energy of the intermediate state as h(E). 

In Fig. [21 we show calculated results for hui\ as a func- 
tion of E F for a (10, 10) armchair nanotube. Here we take 
1595 cm" 1 and 1610 cm" 1 for hLj { x ] of the TO and LO 
modes, respectively. The energy bars denote Y\ values. 
We have used the extended tight-binding scheme to cal- 
culate E e (k), Eh(k), and the electron wavefunction for 
(eh(k)|Wj n t|wA); — As for the el-ph matrix element,— we 
adopted the deformation potential derived on the basis 
of density-functional theory by Porezag et al£L We show 
the resulting Tioj\ as a function of E F at the room tem- 
perature (T = 300 K) and T = 10 K in Fig. [^a) and 
(b), respectively. It is shown that the TO mode does 
not exhibit any energy change while the LO mode shows 
an energy shift and broadening. As we have mentioned 
above, the minimum energy is realized at |-Ef| = hcj^ /2 
(« 0.1 eV). There is a local maximum for the spectral 
peak at I-EfI = 0. The broadening for the LO mode 
occurs within I^f) ^ Huj^ /2 for the lower temperature 
while the broadening has a tail at room temperature for 
\E F \ > Tiuj^/2. 

A continuous model for electrons in a carbon nanotube 
is adopted in this paper to explain the lack of an energy 
shift of the TO modes for armchair nanotubes. As we 
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FIG. 2: (color online) The Ef dependence of the LO (red 
curve) and TO (black curve) phonon energy in the case of the 
(10, 10) armchair nanotube. (a) is taken at room temperature 
and (b) is at 10 K. Only the energy of the LO mode is shifted, 
with the TO mode frequency being independent of Ef. The 
decay width (Fa) is plotted as an error-bar. 



will show in Sec. IIII1 the el-ph matrix element for the 
electron-hole pair creation by the LO and TO modes is 
given by 



(eh(k)|W int |wL > 
(eh(k)|H int | WT o) 



-igu sin#(k), 
-igu cos#(k). 



(2) 



where u is the phonon amplitude, and g is the el-ph cou- 
pling constant. Here 0(k) is the angle for the polar co- 
ordinate around the K (or K') point in the 2D BZ, in 
which a k = (ki, k%) point on a cutting line for a metallic 
energy subband is taken. The k\ (fa) axis is taken in the 
direction of the nanotube circumferential (axis) direction 
(see Fig. |3J|. Equation (JSJ shows that (eh(k)|7Yi n t|w A ) 
depends only on 0(k) but not on |k|, which means that 
the dependence of this matrix element on E is negligible. 
Since the armchair nanotube is free from the curvature 
effect, 22 the cutting line for its metallic energy band lies 
on the &2 axis. Thus, we have 8(k) = ir/2 (— 7r/2) for the 
metallic energy subband which has k\ = and k 2 > 
(k 2 < 0). Then, Eq. © tells us that only the LO mode 
couples to an electron-hole pair and the TO mode is not 
coupled to an electron-hole pair for armchair SWNTs. 




Cutting line 



FIG. 3: (a) Cutting line near the K-point. The ki (k2) axis is 
selected as the nanotube circumferential (axis) direction. The 
amplitude for an electron-hole pair creation depends strongly 
on the relative position of the cutting line from the K-point. 
(b) If the cutting line crosses the K-point, then the angle 
0(k) (= arctan(fc2/fci)) takes n/2 (— 7r/2) values for k-z > 
{ki < 0). In this case, the LO mode strongly couples to an 
electron-hole pair, while the TO mode is decoupled from the 
electron-hole pair according to Eq. ([2]). 



In Fig. 0Ja), we show calculated results for fiuj\ as a 
function of Ep for a (15,0) metallic zigzag nanotube. In 
the case of zigzag nanotubes, not only the LO mode but 
also the TO mode couples with electron- hole pairs. The 
spectrum peak position for the TO mode becomes harder 
for Ep = 0, since Ke(h(E)) for E < S-wto contributes 
to a positive frequency shift. It has been shown theo- 
retically 22 and experimentally 23 that even for "metallic" 
zigzag nanotubes a finite curvature opens a small energy 
gap. When the curvature effect is taken into account, the 
cutting line does not lie on the K-point, but is shifted 
from the ki axis. In this case, cos#(k) = k\j(k\ -Ffc 2 ) 1 / 2 
is nonzero for the lower energy intermediate electron-hole 
pair states due to k\ ^ 0. Thus, the TO mode can cou- 
ple to the low energy electron-hole pair which makes a 
positive energy contribution to the phonon energy shift. 
The high energy electron-hole pair still decouples to the 
TO mode since cos6*(k) — * for \k%\ ^> Therefore, 
when \E-p\ ^5 then Ti^to increases by a larger 

amount than Tiujlo- The TO mode for the small diam- 
eter zigzag nanotubes couples strongly with an electron- 
hole pair because of the stronger curvature effect. In 
Fig. 0Jb), we show the diameter (dt) dependence of the 
fiij\ of zigzag nanotubes for Ep = not only for metallic 
SWNTs but also for semiconducting SWNTs. In the case 
of the semiconducting nanotubes, the LO (TO) mode 
appears around 1600 (1560) cm -1 without any broaden- 
ing. Only the metallic zigzag nanotubes show an energy 
shift, and the energy of the LO (TO) mode decreases 
(increases) as compared to the semiconducting tubes. In 
the lower part of Fig. [4^b), we show a curvature-induced 
energy gap E gap as a function of dt- The results show 
that higher (lower) energy electron-hole pairs contribute 
effectively to the LO (TO) mode softening (hardening) in 
metallic nanotubes. In the case of semiconducting nan- 
otubes, we may expect that there is a softening for the 
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LO and TO modes according to Eq. ([2]). However, the 
softening is small as compared with that of the metallic 
nanotubes because the energy of intermediate electron- 
hole pair states is much larger than in this case. 
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FIG. 4: (color online) (a) The Ef dependence of the LO (red 
curve) and TO (black curve) phonon frequency for a (15, 0) 
zigzag nanotube. Not only the frequency of the LO mode but 
also that of the TO mode is shifted due to the curvature ef- 
fect, (b) The diameter dt dependence of the phonon frequency 
for zigzag nanotubes, including zigzag semiconducting tubes. 
.Egap denotes the curvature-induced mini energy gap. 



electrons with energy eigenvalue E is written by 
v f ct ■ (p + A(r))#(r) = £tf (r), 



(3) 



where vf (= Jo^/K) is the Fermi velocity and 70 (~ 
2.7eV) is the nearest-neighbor hopping integral, £ = 
(3/2)a cc , p = —ihSJ is the momentum operator, and 
er = (<JxjOy) is the Pauli matrix. A(r) is given in Eq. (3) 
of Ref. 25 by a small change S^f a (r) (a = 1,2,3) of the 
hopping integral from —70 (see Fig.[S]), as 



v F A x (r) = &yfo) - \ (<5 7o 2 (r) + S^(t)) 
^A,(r) = ^(,57o 2 (r)-57o 3 (r)). 



(4) 



Here ^7g( r ) for the LO and TO modes is given by 
<$7o( r ) = (ff/^) u ( r ) • R-a where R a denotes the nearest- 
neighbor vectors (Fig. [5|) and u(r) is the relative dis- 
placement vector of a B site from an A site (u(r) — 
Us(r) — Ujv(r)) and g is the el-ph coupling constant. We 
rewrite Eq. Q as 



v F (A x (r),Ay(r)) = g(u y (r),-u x (r)), 



(5) 



x,y), and Ri-(R 2 -rR 3 )/2 = 
— £e x have been used (see 



where iii(r) = u(r)-e.;, (i 
le y and V3/2(R 2 - R3) = 
Fig. [5]) . Then, the el-ph interaction for an in-plane lattice 
distortion u(r) can be rewritten as the vector product of 
er and u(r))i£ 



vfct ■ A(r) = g(cr x u(r)) ■ e 2 



(6) 




III. CHIRALITY DEPENDENCE OF THE 
ELECTRON-PHONON INTERACTION 

Here, we derive Eq. ([2|) on the basis of a continu- 
ous model for 7r-electrons near the K-point in graphcnc. 
In a continuous model, the local modulation of the 
hopping integral due to lattice vibrations appears as a 
deformation-induced gauge field, A(r) = (A x (r),A y (r)), 
in the Weyl equation^ The Weyl equation for ir- 



FIG. 5: A hexagonal unit cell of graphene consists of A (closed 
circle) and B (open circle) sublattices. R a (a = 1,2,3) are 
vectors pointing to the nearest-neighbor B sites from an A site 
(Ri = a cc e y , R2 = —(V3/2)a cc e x — (l/2)a cc e y , and R3 = 
(V3/2)a cc e x — (l/2)a cc e v ). Local modulations of the hopping 
integral are defined by #70 (r) (a = 1, 2, 3). The modulation is 
given by optical phonon modes as £70 = (g/£)u(r) ■ R a where 
u(r) (= ua(i") — ua(r)) is a relative displacement vector of a 
B site relative to the nearest A site. 

We consider the LO and TO phonon modes with q = 
(i.e., T-point). Then, an electron-hole pair is excited by 
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a constant u = (u x ,u y ). The el-ph matrix element for 
the electron-hole pair generation is given by 

(eh(k)\n int \u) = J *; k (r)W int ^ ik (r)d 2 r 

/ \ / - e(k) \ 

5 / 4., e < k ) ,-Q( k A / Uy + iu x \ e % 2 \ 
= 2 l e+ 2 6 2 J \u y -iu x J \_ e -H^J 

(7) 

where vE , c ,k(r) (^«,k(r)) denotes an energy eigenstate of 
vf<j ■ p in the conduction (valence) energy band with 
energy eigenvalue E — «f|p| (E = — Vp\p\), 

gik-r / -ie(k)/2\ 

e ik-r / e -i0(k)/2 \ ^ > 

^ k(r) = 7H U +i6(k)/ V' 

where 5 denotes the total area of the system, k x — ik y = 
|k|e- ie ( k ). 

We first consider the case of a zigzag nanotube in 
Fig. [5] Then, we denote x (y) as a coordinate around 
(along) the axis (so that 6(k) = 0(k)), and ^(r) (u y (r)) 
are assigned to the TO (LO) phonon mode. The corre- 
sponding A(r) for u x (r) and u y (r) is A y (r) and A x (r), re- 
spectively. By calculating Eq. ([7]) for the TO mode with 
(u x ,u y ) = (u,0) and for the LO mode with (u x ,u y ) — 
(0,u), we get Eq. ©. Next, we consider the case of 
an armchair nanotube. Then, x (y) is the coordinate 
along (around) the axis (so that 0(k) = 6(k) +n/2), and 
u x (r) (u y (r)) is assigned to the LO (TO) phonon mode. 
The direction of the gauge field A(r) is perpendicular to 
the phonon eigenvector u(r) and the LO mode shifts the 
wavevector around the tube axis, which explains how the 
LO mode may induce a dynamical energy band-gap in 
metallic nanotubes By calculating Eq. ([7]) for the TO 
mode with (u x ,u y ) — (0,it) and for the LO mode with 
(u x ,u y ) = (u, 0), we get Eq. ((2|), too. Equation ((2|) is 
valid regardless of the tube chirality if the phonon eigen- 
vector of the LO (TO) phonon mode is in the direction 
along (around) the tube axis. This is because p and u(r) 
are transformed in the same way as we change the chi- 
ral angle. As a result, there would be no chiral angle 
dependence for the el-ph coupling in Eq. ((2]). 

However, the phonon eigenvector depends on the chiral 
angle. Reich et al. reported that, for a chiral nanotube, 
atoms vibrate along the direction of the carbon-carbon 
bonds and not along the axis or the circumference pi£ In 
the case of a chiral nanotube, the phonon eigenvector 
may be written as 

fu TO \ _( cos0 sin (A /VA , , 

V"lo/ \-sin0 cos 4>) \U2j' 

where u\ (1x2) is in the direction around (along) a chiral 
tube axis, and 4> is the angle difference between the axis 
and the vibration. This modifies Eq. @ as 

(eh(k)|Hint|wLo) = -igusm(9(k) + <f>), 
(eh(k)|Hi n tkTo) = -igu cos(#(k) + 0). ' 



In Fig. [6l we show numerical results for T%uj\ as a function 
of Ep for a (10, 4) metallic chiral nanotube for T = 300 
K. The energy difference between the minimum point at 
Ey = 7kj'°) /2 and at E-p = shows that the LO mode 
couples more strongly to the low energy electron-hole pair 
than the TO mode. Since the curvature-induced mini en- 
ergy gap for (10, 4) and (15, 0) tubes are almost the same 
(.©gap = 0.045 eV), we may expect a similar energy shift 
of the LO and TO modes. However, the results of the 
(15, 0) zigzag tube show that the TO mode couples more 
strongly for the low energy electron-hole pair than the LO 
mode. This is explained not by Eq. @ but by Eq. fTU]) 
with appropriate (j> ((f) 23.4°), which is given by the 
phonon eigenvector calculation. Since the chiral angle 
for (10,4) is 16.1°, <j> is not directly related to the chiral 
angle. The identification of cf> in Eq. (fTT)]) as a function of 
chirality would be useful to compare theoretical results 
and experiments, which will be explored in the future. 



T = 300 K £ gap = 0.045 eV 
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FIG. 6: (color online) The Ef dependence of the LO (red 
curve) and TO (black curve) phonon energy for the (10, 4) 
metallic chiral nanotube. -E gap denotes the curvature-induced 
mini energy gap which has the same value as that for the 
(15, 0) tube in Fig.|4|a). The difference between the behavior 
of the (10, 4) and (15, 0) nanotubes comes from the fact that 
the phonon eigenvector of the LO (TO) mode is not pointing 
along the tube axis in the case of the (10, 4) tube. 



IV. DISCUSSION 

The phonon frequencies with some broadening as 
shown in Fig. [2] are not directly related to the Raman 
spectra of the G-band but relate to the phonon density 
of states of the LO and TO phonons at q = 0. In previ- 
ous papers^ we have shown that the G-band intensity 
depends on chiral angle in which for zigzag nanotubes 
the LO (TO) phonon mode is strong (suppressed) while 
for armchair nanotubes the TO (LO) phonon mode is 
strong (comparable). This chirality dependence of the 
G-band intensity is also observed in single nanotube Ra- 
man spectroscopy ! 27 ! 28 This chirality dependence of the 
G-band mode intensity is exactly opposite to the chi- 
rality dependence of the phonon-softening in which the 
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TO phonon mode is suppressed in armchair nanotubes. 
These observations clearly show that the corresponding 
el-ph interaction for phonon-softening and Raman pro- 
cesses are independent of each other. 

In the Raman process, the phonon is emitted by scat- 
tering an electron (or an exciton) from one conduction 
band state to another conduction band state, while in 
the phonon softening, the electron is scattered (or ex- 
cited) from a valence band state to a conduction band 
state. Thus the corresponding matrix element for the 
Raman process is expressed by substituting x P Cj k( r ) m 
Eq. |n|) for *«,k(r) in Eq. fl7j as follows, 

(e(k),w|W int |e(k)) = J * c * k (r)H in t* c ,k(r)cPr 

9 ( , .eoo _,e<k)\ / u y + tu x \ / e 1 2 \ 
= 2 l e+ 2 6 2 ) [uy-iu x ) [ e+ i^)- 

(11) 

Then the electron-phonon matrix element for the Raman 
scattering process becomes 

(e(k),w LO |Hint|e(k)) = .gucos6»(k), 
(c(k),c^ro|Wint|e(k)> = ~gusm9{k). 

It is stressed again that the el-ph interactions of Eqs. ((2|) 
and (|12p are for phonon-softening and Raman inten- 
sity, respectively and that they are independent of each 
other. When we consider the Raman scattering processes 
for carbon nanotubes, the k vector which is relevant to 
the resonance Raman process is the k vector at the van 
Hove singular point, k&, which corresponds to the touch- 
ing points of the equi-energy surface to the cutting line 
(one dimensional Brillouin zone plotted in the 2D BZ)i^ 
which are shown in Fig. [71(a) and (b) for zigzag and arm- 
chair nanotubes, respectively. 




Armchair (9 fO) Zigzag (6 = 0) 



FIG. 7: Trigonal warping effect of an equi-energy surface in 
the 2D Brillouin zone and cutting lines for (a) armchair and 
(b) zigzag nanotubes. In the case of (a) armchair nanotubes, 
the 9 value defined from the fci axis in Fig. [5] is not zero 
while in the case of (b) zigzag nanotubes, 9 is always zero, 
which explains the chirality dependence of the relative Raman 
intensity of LO and TO phonon modes. 

In the case of the zigzag nanotubes, the k,j point cor- 
responds to 9 = which is not the chiral angle, but is 
defined in Fig. [3](b). Thus from Eq. (|12p. we see that 



the TO mode is not excited. On the other hand, in the 
case of the armchair nanotube, the value of 9 is not zero 
because of the trigonally warped equi-energy surface as 
shown in Fig.[7£b)iii22 Thus not only the LO but also the 
TO phonon modes can be observed, which is consistent 
with the previous experimental result o 27 i 28 and the the- 
ories'^ Thus when we see the G-band Raman spectra 
for metallic nanotubes, we can expect the following fea- 
tures: for armchair nanotubes, we can see both LO and 
TO in which the LO mode appears at a lower frequency 
region with some spectral width while the TO mode ap- 
pears without spectral broadening because of the absence 
of phonon softening. For zigzag nanotubes, we can see 
only a broadened LO phonon mode but we can not see 
a TO phonon mode. For a general chiral nanotube, we 
generally observe both LO and TO phonon modes with 
some broadening as a function of 0. Further, we expect 
some phonon hardening effect for the TO phonon mode. 
The relative intensity between the LO and TO phonon 
modes is determined by the chirality dependent Raman 
intensity and spectral width. The detailed Raman spec- 
tral features for all metallic (n, m) nanotubes will be pre- 
sented elsewhere. 

The gauge field descriptions for the lattice deformation 
(Eq. Q) and for the el-ph interaction of the LO and TO 
modes (Eq. |5])) are useful to show the appearance of 
the curvature-induced mini energy gap in metallic zigzag 
carbon nanotubes and the decoupling between the TO 
mode with a finite wavevector and the electrons, as shown 
in the following. For a zigzag nanotube, we have S^/q = 
and ^7g = Sjq from the rotational symmetry around 
the tube axis (see Fig. [5}. Then, Eq. ([4]) shows that for 
A x ^ and A y = 0, the cutting line of k x = for the 
metallic zigzag nanotube is shifted by a finite constant 
value of A x because of the Aharanov-Bohm effect for the 
lattice distortion-induced gauge field A. This explains 
the appearance of the curvature- induced mini energy gap 
in metallic zigzag carbon nanotubes^ and of the phonon 
broadening for the TO mode as a function of E-p. The 
TO phonon mode with q / does not change the area 
of the hexagonal lattice but instead gives rise to a shear 
deformation. Thus, the TO mode (uxo( r )) satisfies 

V • u TO (r) =0, V x u TO (r) + 0. (13) 

Using Eqs. ([5]) and (fTB")) . we see that the TO mode does 
not yield the deformation-induced magnetic field (B(r) = 
V x A(r)), but the divergence of A(r) instead because 

B z (r) = — ^V-u T o(r) = 0, 

(14) 

V ■ A(r) - ^- V x Uto r) ■ e z ^ 0. 

Thus, we can set A(r) = Vy>(r) where ip(r) is a scalar 
function. Since we can set V • A(r) = in Eq. © by 
a redefinition of the phase of the wavefunction (a local 
gauge symmetry) as ^(r) — ► exp(— iip(r)/h)^(r)^ and 
thus the A(r) in Eq. ([3]) disappears for the TO mode 
with q^0. This explains why the TO mode with q ^ 
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completely decouples from the electrons and that only 
the TO mode with q = couples with electrons. This 
conclusion is valid even when the graphene sheet has a 
static surface deformation. In this sense, the TO phonon 
mode at the T-point is anomalous since the el-ph interac- 
tion for the TO mode can not be eliminated by a phase of 
the wavefunction. It may be difficult to include the local 
gauge symmetry in a numerical analysis. On the other 
hand, the LO phonon mode with q ^ changes the area 
of the hexagonal lattice while it does not give rise to a 
shear deformation. Thus, the LO mode (uLo( r )) satisfies 



V-u LO (r)^0, Vxu LO (r)-0, 



(15) 



Using Eqs. ([5]) and (fT5"j) . we see that the LO mode gives 
rise to a deformation-induced magnetic field as 



B z (t)^0, V-A(r)=0. 



(16) 



Since a magnetic field changes the energy band structure 
of electrons, the LO mode (with q ^ 0) couples strongly 
to the electrons. 

In the case of 2D graphene, Eq. ([2]) tells us that the 
L point TO and LO modes give the same energy shift 

(2) 

because the integral over #(k) gives the same ?kj^ in 
Eq. ([1} for both TO and LO modes. This explains 
why no G-band splitting is observed in a single layer 
of graphene^i Even when we consider the TO and LO 
modes near the L point, we do not expect any splitting 
between the LO and TO phonon energies because the TO 
mode with q ^ is completely decoupled from the elec- 
trons. Thus, for q ^ 0, only the LO mode contributes to 
the G-band. 

In our numerical results for Tiuj\, we see a crossing be- 
tween the LO and TO mode at some values of E-p. If 
there is a finite transition amplitude between the LO and 
TO phonon modes, we need to diagonalize the following 
matrix according to second-order perturbation theory for 
a degenerate case, 

/>L |H (2) | 



|wto ; 



where HS 1 ^ is the effective Hamiltonian in second-order 
perturbation theory given by 



E 



H in t|eh(k))(eh(k)|H in t 



fiuf - (E c (k) - £ h (k)) + if 



(18) 



Using Eq. ([2]), we see that cos#(k) sin6'(k) which is an 
odd function of ki appears in the numerator for the diag- 
onal terms in Eq. (|18p . Thus, the mode coupling through 
the electron-hole pair at ki is canceled by that at — 
i.e., due to the time-reversal symmetry of the system. 
Thus, even though the LO and TO modes cross each 
other, there is no-mixing between the LO and TO phonon 
modes for any chirality. If there is a breaking of time re- 
versal symmetry for fc 2 , we expect some mixing between 
the LO and TO phonon modes. 



V. COMPARISON WITH EXPERIMENT AND 
SUMMARY 

The first comparison we make is to experimental re- 
sults in the literature. Here we try to assign the chiral- 
ity of the nanotube that is given in Fig. 2 in Ref. HtH 
This figure shows two strong intensity peaks with differ- 
ent phonon frequencies. The higher frequency peak does 
not depend on the gate voltage and the lower frequency 
peak shows a frequency shift with a strong broadening 
near the Dirac point. The existence of a flat intensity 
peak as a function of Fermi energy for a metallic SWNT 
indicates that the electron-phonon coupling is very weak 
for the phonon mode. Thus, the nanotube can be thought 
of as an armchair SWNT (or close to an armchair SWNT 
in the chiral angle) since the el-ph coupling between the 
TO mode and electron-hole pairs is negligible in this case. 
Further investigation comparing theory with experiment 
is strongly necessary for tubes of different chiralities. 

In Fig. [5Ja) and (b), we show experimental results of 
the G-band intensity as a function of applied gate voltage 
for two different isolated SWNTs. For the given excita- 
tion energy of 1.91eV and the observed radial breathing 
mode (RBM) frequencies 196 and 193 cm -1 for the sam- 
ple of Fig.[8]Ja) and (b), respectively, we can assign (n, m) 
by using the conventional assignment technique that we 
used for single nanotube Raman spectroscopy. 32 As a re- 
sult, we assign the (n, m) values as (12,6) for Fig. [5] (a) 
and (15,3), (16,1), or (11,8) for Fig.^b). 



(a) 



1650| 



2 
Gate Voltage [VI 




-1 -0.5 0.5 
Gate Voltage [V] 



FIG. 8: (color online) Experimental results of the Raman G- 
band intensity as a function of applied gate voltage for two 
metallic SWNTs (a) and (b). A strong (weak) intensity peak 
is denoted by the yellow (black) color. 

Although (12,6) is not so close to the chiral angle for 
the armchair nanotube, we expect that Figure [HI a) ex- 
hibits a similar behavior to that of Fig. 2 in Ref.[l6|. Thus 
the assignment of (12,6) from the RBM spectra is consis- 
tent with the present phonon spectra. In Fig. [HJb), the 
intensity near the Dirac point of the TO (LO) mode is 
weaker (stronger) than shown in (a). This indicates that 
the tube is not an armchair SWNT but close to a metal- 
lic zigzag SWNT (or a metallic chiral SWNT). Thus we 
can assign this SWNT either to (15,3) or to (16,1) for 
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the SWNT for Fig. [gjb) and we can exclude (11,8). A 
further comparison with more experimental data will be 
required, which will be carried out in the future. 

A chirality dependent Raman intensity was observed 
for metallic SWNTs as shown at Fig. 1 in Ref. [H. 
Although the Ep positions for the observed metallic 
SWNTs are unclear, the results are consistent with our 
calculations in the following way. For example, the TO 
mode in a (15, 15) SWNT gives a sharp lineshape and the 
LO mode shows a broad feature for a (24, 0) SWNT. It is 
noted that the Raman intensity is proportional to the el- 
ph coupling for an optically excited electron via Eq. (fT2")l . 
Since this el-ph coupling depends on the chiral angle due 
to the trigonal wrapping effect, one of the two optical 
modes may be invisible in the Raman intensit y 33 ' 34 

In summary, we have calculated the E-p dependence of 
hcu\ (A = LO and TO) for metallic carbon nanotubes. 
The results show a strong dependence of the phonon 
frequency shift on the chirality of single-walled carbon 
nanotubes because of the curvature-induced shift of the 



wavevector around the tube axis. This is explained by 
the general property of second-order perturbation theory 
and the characteristic electron-phonon coupling for the 
T point LO and TO phonon modes (Eqs. and ©). 
For the LO and TO phonon modes near the T point, 
we showed the absence of an electron-phonon coupling 
for the TO mode for q ^ due to a local gauge sym- 
metry, and that the LO mode works as a deformation- 
induced magnetic field. The phonon mixing between LO 
and TO phonon modes is absent in second-order pertur- 
bation theory due to time-reversal symmetry. 
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